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' Primordial gravitational waves are created during the de Sitter phase of an 

^ ■ exponentially-expanding (inflationary) universe, due to quantum zero-point 

. . vacuum fluctuations. These waves produce fluctuations in the temperature of 

the Cosmic Background Radiation (CBR). We calculate the multipole mo- 
ments of the correlation function for these temperature fluctuations in a 
spatially-closed Friedman-Robertson- Walker (FRW) cosmological model. The 
results are compared to the corresponding multipoles in the spatially-flat case. 
The differences are small unless the density parameter today, Qq, is greater 
than 2. 

PACS numbers: 98.80.Cq, 98.80.C, 98.80.Es 
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I. INTRODUCTION 



Inflationary models of the early universe contain a well-studied mechanism which creates 
primordial fluctuations. The fluctuations originate as quantum-mechanical zero-point fluc- 
tuations during the exponentially-expanding de Sitter phase. By a process which may be 
variously described as particle (graviton) production, non-adiabatic ampliflcation, or super- 
radiant scattering, these fluctuations become large in the present epoch. As the universe 
expands, these perturbations are redshifted to longer wavelengths and amplifled; during 
the present epoch these perturbations typically persist over a range of wavelengths A from 
10^^'' cm < A < 10^ cm. For a review of perturbations in inflationary models, see Kolb and 
Turner 

The perturbations of the gravitational fleld may be decomposed into scalar, vector and 
tensor components. The tensor perturbations considered in this paper may be thought of 
as gravitational waves in a classical description, or as spin-two gravitons in the quantum 
mechanical description used in the present work. The modes of interest have present-day 
frequencies in the range from 10^^^ Hz to 10~^^ Hz and have extremely large occupation 
numbers. Hence they may also be thought of as classical gravitational waves - the two 
descriptions coincide. The gravitons are created during the de Sitter phase of rapid expansion 
by the mechanism originally proposed by Parker; the same mechanism creates particles near 
a black hole or in any other region where the spacetime curvature is large and particle 
creation is not forbidden by global symmetries or conservation laws. A simple calculation 
showing how a potentially-observable spectrum of gravitons is created in inflation is given 



The tensor perturbations of the gravitational fleld produce temperature fluctuations in 
the CBR, via the Sachs- Wolfe effect. The expected values of the resulting temperature 
fluctuations are described by the angular correlation function 



Here 6T/T{Q) is the fractional temperature fluctuation in the CBR at point Q on the 
observer's celestial sphere, 7 is the angle between Q and Q', and the quantum expectation 
value is evaluated in the initial state of the universe. The multipole moments {af) are 
generally used to describe C(7). 

In a recent paper the expected multipole moments {af) due to tensor perturbations 
are calculated in a spatially flat k = FRW inflationary model. That paper contains a 
detailed review of previous work on this problem, a comprehensive description of the physical 
motivation, and a detailed and self-contained "flrst-principles" calculation. The present work 
repeats that calculation in the spatially-closed {k = +1) case. The only previous work on 
tensor perturbations in the spatially-closed case is that of Abbott and Schaefer [Q. Note 
that the angle brackets around af serve as a reminder that we are calculating the expected 
or expectation values of these multipole moments, not necessarily the values that they might 
have in any given realization of the universe. 

The calculation in this paper follows the previous work by Allen and Koranda ^ very 
closely. In the present work, we will assume that the reader is familiar with that earlier 
paper, and present only the bare minimum of detail required to generalize the work to 
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the k = +1 case. In Section |T| we present the k = +1 cosmological model and Sachs- 
Wolfe effect. Section III gives the form of the metric perturbation operator for linearized 
gravitational fluctuations. Section |IV| combines these results to obtain an analytic form 
for the multipole moments (af). Section ^ details the method by which these multipole 
moments were evaluated numerically, and Section |VT| outlines the results and conclusions of 
that numerical study. 

Throughout this paper, we use units where the speed of light c = 1. However for clarity 
we have retained Newton's gravitational constant G and Planck's constant h explicitly. We 
choose function branches so that y/x > and arcsin(x) G [— 7r/2, 7r/2]. 



II. THE BACKGROUND SPACE-TIME AND THE SACHS- WOLFE EFFECT 

The spacetime considered here has the topology R x of the static Einstein cylinder, 
and is covered by coordinates = t,x^ = x^^"^ = ^^^^ = with the ranges X; ^ ^ [0,7r], 
and G [0, 2n). The time coordinate t ranges over a connected open subset of the real line, 
which we will specify below. The spatial coordinates cover a three-sphere of radius a{t); we 
refer to this function as the cosmological scale factor. The metric of the spacetime is given 
by 



ds^ = a\t) l^-dr + dx + sin^ xidO^ + sin^ dd(f)^) + Kjit, x, ^, <p)dx'dx^y (2.1) 

The metric perturbation hij is assumed to be small; in its absence the spacetime metric is 
that of a homogeneous and isotropic k = +1 FRW model. With our choice of gauge for the 
tensor metric perturbations, the indices i,j = 1,2,3 run only over the spatial coordinates. 

In order to completely specify the cosmological model, we need to define the cosmological 
scale-factor a{t). The cosmological model is completely defined by the free parameters given 
in Table |. Note that we have assumed that the universe is currently expanding, since we 
require Hq to be positive. The density parameter 

is the ratio of the present-day energy-density po to the critical energy density required to 
produce a spatially-fiat k = universe. 



A. The Matter-Dominated (Dust) Phase 

In our cosmological model, the universe is assumed to pass through three "phases", 
appropriate to a simple inflationary model. We let t = denote the present time. The most 
recent phase was a matter-dominated period of expansion, described by the scale factor 

a{t) = Asm^{t/2 + B) for teq < ^ < 0. (2.3) 

Here the constants A, B, teq are defined by 
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-no 



i? = arcsin W — — , and (2.4) 

V "0 



teq = 2 arcsm ^ / ^^^^ ^r—. — 2 arcsm 



noil + Zeq) V ^0 



During this matter-dominated phase, the stress-energy tensor is that of a perfect fluid, with 
zero pressure and an energy density proportional to a~^{t). We assume (as indicated in 
Table |) that the surface of last scattering is located within the matter-dominated phase. 
Thus, the time of last scattering. 



is given by a formula identical in form to ( ^.41 ) for teq, and satisfies teq < tis < 0. 

B. The Radiation-Dominated Phase 

Preceding the matter-dominated phase of expansion is a radiation-dominated phase of 
expansion. During this phase the scale factor is 

a{t) =Csm{t + D) for tend < t < teq. (2.6) 

Here the constants C, D, tend are defined by 



rin 



iJ = 2arcsmw — — arcsm w-— — — , and 2.7 

V V^o(l + ^cq) 



fend = — 2 arcsm \l — — h arcsm ^ / ^ — ■ — ^^-^ -|- arcsm , 



no V^0(l + ^eq) \ ^oil + Z, 



(^]o-l)(l + ^, 



cq; 



^cnd)^ 



During this radiation-dominated phase of expansion the energy density is proportional to 
a~'^(t) and the pressure is equal to 1/3 of the energy-density. This phase is preceded by a 
de Sitter phase. 

C. The Initial de Sitter (Inflationary) Phase 

In our coordinate system, the de Sitter (exponentially expanding, inflationary) phase has 
scale factor 

a{t) = for t<tend. (2.8) 

sm(t + r ) 

Here the constants E and F are defined by 
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E = + Z,^Y/\l + Zend)-', and (2.9) 

-no 



r = 2 arcsm \ / — — arcsm \ / 7^— --— — 2 arcsm 



fio V^0(l + ^eq) \ ^^0(1 + ^, 



(^]0-l)(l + ^eq) 



^end)^ 



Note that the constant E < because sin(t + F) < during the de Sitter phase. During 
the de Sitter phase, the energy density is a constant 

_ 3 2 _ ^ (1 + Zcnd)^ _ 3i/o^o (1 + ^cnd)^ .^v 
PdeSittcr - — Po -I , ^ - g rT~V l^'-l-Uj 

and the (negative) pressure is — pdeSitter- 

D. Properties of the Cosmological Model 

It may be easily verified that the scale factor and its derivative w.r.t. time t are both 
continuous, however the second derivative is discontinuous. This is because in our simple 
infiationary model, the energy-density is a continuous function but the pressure changes 
discontinuously at the beginning and end of the radiation-dominated epoch. 

The de Sitter phase "begins" at early times when the time coordinate t approaches the 
value —71 — F. At this early time the cosmological scale factor is very large (approaching 
infinity as t ^ — vr — F). As the time coordinate increases, the scale factor decreases, 
eventually reaching a minimum value when t = tmin = — 7r/2 — F. After this time, the 
scale factor begins to increases again (exponentially in physical time). One might find it 
reasonable to demand that the universe be expanding at time tend when the infiationary 
phase ends. This is the case if and only if tmm < tend, which implies that the free parameters 
given in Table |I| must satisfy the inequality 



/(fi0-l)(l + Zeq) 



<l + ^end. (2.11) 



It is also easy to determine the "amount" of infiation that takes place. The amount that 
the universe has expanded between time tmin, when the spatial sections have their smallest 
extent, and time tend, when the infiationary phase terminates and the radiation-dominated 
phase begins, is 



o(tcnd) n , ry \ ^0 



a(tmin) ' ^'="""V(^0-1)(1 + ^, 



(2.12) 



Comparison with (|2.11| ) shows the obvious - if the universe is expanding at the end of the 
de Sitter phase, then the amount of infiationary expansion (|2.12| ) is greater than unity. In 
typical infiationary models, the free parameters have values of order Hq between 50 and 100 
Km/s-Mpc, fio < 2, 100 < Zi, < 1500, 2 x 10^ < Zeq < 2 x 10^ and lO^^ < Zend- 
There is a sense in which the spatially-closed infiationary models are not "natural." One 
of the principal motivations which led to the development of the infiationary paradigm was 
the desire to solve the so-called "horizon problem." As we will now show, this problem is 
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only solved (for reasonable choices of the cosmological parameters) if < 2. Thus, while 
it is technically consistent to use the results obtained in this paper for any value of fio > 1, 
one must bear in mind that the cosmological model, for large values of Qq, runs counter to 
the spirit of inflation. 

The horizon problem may be stated in terms of a set of points C, which is the intersection 
of the past horizon of an observer today with the surface of last scattering. The horizon 
problem is "solved" if C lies within the causal domain of influence of either (1) a point on 
the initial singularity, in a big bang model, or (2) a point on the surface at t = t^i^ where 
inflation "begins", in a model with no initial singularity. Thus, in our model, which is of 
type (2), the horizon problem is solved if and only if 



1^0 — ^ls| < l^ls — t 



mm 



(2.13) 



2 arcsin , - 4 arcsin , /^^Wf^ + arcsin , f^W^f^ + 2 arcsin , h^^P.'-l 



a ^ ± 

no V f^0(l + ^ls) ' """""" V f^0(l + 2cq) ' Y no(l + Zend)2 2' 



For reasonable cosmological models, the terms containing Zi^, Z^q, and Z^nd niay be ne- 
glected. The horizon problem is then solved if and only if 




arcsin W^^- — < - fio < 2. (2.14) 



While we present results for any value of Qq, the cosmological model itself should be viewed 
with some suspicion if Qq is much larger than unity. 

The final result of this paper are values of the dimensionless quantities 

Ml ^ ^"^^^ + ^\ af). (2.15) 

PdeSittcr 6 

Here ppianck is the Planck energy-density ppianck = ~ 5 x 10^'^ gm/cm^. It will turn out 
that Ml is independent of Hq, and depends only upon the dimensionless quantities Qq, Z\s, 
Zeq, and Zend- This is because pdeSitter, and (of) are both proportional to Hq. In addition, 
if Zend is sufficiently large then the Mi are also independent of its value. 



E. The Sachs- Wolfe EfTect 

If the metric perturbation hij vanishes and the temperature of the CBR on the surface 
of last scattering is constant, an observer today would see exactly the same temperature at 
each point on the celestial sphere, and C(7) would vanish. However the metric perturbations 
will in general break the rotational symmetry and perturb the energy of the photons. This 
results in a temperature fluctuation which varies from point to point on the celestial sphere; 
the fluctuation may be calculated in the same way as for a spatially-flat Universe, given in 

i- 

We assume that the observer is located at t = 0, and at "radial" coordinate x = 0. 
(Because the coordinate system is singular at x = every value of 6, corresponds to the 
same space-time point at x = 0, so their values are irrelevant when x = 0-) If the observer 
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looks out at a point Q on the celestial sphere, she observes photons that arrive, in the 
unperturbed metric, along the null geodesic path 

t(A) = A x(A) = |A| = -A e{\)=9n = (pn- (2.16) 

In these equations, 6^ and (pn are the angular coordinates of the point Q on the celestial 
two-sphere. We have chosen the (non-affine) parameter A along the null geodesic path to 
run through the range tis < A < between the time of last scattering and the observation 
today. 

In the presence of the metric perturbation hij the fractional temperature fluctuation 
observed at point Q on the celestial sphere is 

ST 1 /-o /dh \ 

= -2 dX j (t(A), x(A), ^(A), 0(A)). (2.17) 

As indicated in this formula, the partial derivative w.r.t. the time coordinate t is taken 
before setting the coordinates equal to the values which they take along the unperturbed 
null geodesic path. 



III. THE METRIC PERTURBATION OPERATOR 



The classical metric perturbation hij may be replaced with a quantum field operator. 
The justification for this is given in our detailed paper on the spatially-flat case and 
will not be repeated here. The basic idea is that the inflationary epoch redshifts away all 
the perturbations, with the exception of the zero-point quantum fluctuations. Hence we 
calculate the expectation value of ^(n)^(n') in the vacuum state |0) appropriate to the 
initial de Sitter state. 



A. Mode function expansion of Metric Perturbation Operator 

The quantum field operator (which we denote with the same symbol h^j as the cor- 
responding classical perturbation) may be expanded in terms of a complete set of mode 
functions. As was originally shown by Ford and Parker 0, in an FRW cosmological model, 
the time-dependence of these mode functions is the same as that of a massless minimally- 
coupled scalar field. The field operator is 



h,,{t,x,OA)= (3.1) 

OD L I 

EE E 



L=2 1=2 m= 



In this expression, the sum is over a complete set of rank-two symmetric transverse traceless 
tensors These tensor modes are defined on a unit-radius sphere 5*^ and are given 

explicitly by Higuchi 0. (Note however a typo [§ in one of the formulae which does not 
affect the results which we need.) Henceforth we will denote the triple sum that appears in 
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by J^Lim without explicitly indicating the ranges of summation. (The temperature T 
can always be distinguished from the tensor modes Tj^^^^\ since the latter is always written 
with indices.) 

The graviton has two possible polarization states, labeled s and v in this expansion, each 
of which has its own set of tensor modes. The modes are labeled by the three integers L, 
I, and m. (Note that Ford and Parker's index n = L + 1 in our notation.) Associated 
with the s-polarization modes are creation and annihilation operators cum and c^um, and 
associated with the f -polarization modes are creation and annihilation operators dum and 
d^;^. The only non-vanishing commutation relation among this infinite set of operators is 
the relation 

[cUm, C^L'l'm'] — [dum, d^L'l'm'] — ^LU^W^mm' (3.2) 

where 5 denotes the Kronecker delta function. 



B. The Transverse- Traceless-Symmetric Tensor Harmonics 

The tensor modes defined by Higuchi [0 obey the normalization condition 

r dxsin^X r define rf^T.^'^^'^^T^P '^'''"^V^^P^'^ = 5LL'5u'6mm'5yp> . (3.3) 

JQ Jo Jo 



Here, the polarization indices p and p' take on either of the values "s" or "v" . The integral 
is over the unit-radius three-sphere, and P^^ is the inverse of the metric on the unit-radius 
three-sphere: 

Pijdx'dx^ = dx^ + sin^ x{dO'^ + sin^ Odcj)'^). (3.4) 
Note that the measure that appears in the normalization integral ( p.3|) is the usual volume 



element defined by ydetPj^. The Sachs- Wolfe effect ( p.l7|) is produced only by the XX 
component of hij. Because T^^'"^'™-) = Q, only the "s" polarization state contributes to the 
temperature fluctuation. The only component needed is thus 

T(f' ™)(x,^,0) = R[{x)YUO,<P). (3.5) 

The Yim{0,(j)) are standard scalar spherical harmonic functions on the two-sphere 0. The 
"radial" dependence is given by 



(/-l)/(/ + l)(/ + 2)(L + l)(L + / + l)! p-(m/2). ^ . .3 

\ 2L{L + lf{L + 2){L-l)\ ^+1/2 icosxj, Ici.bj 



where the functions Pl+i/V^"* (-2) are associated Legendre functions 0. In Section we 
explain how these functions may be easily evaluated. 
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C. Normalization Condition for Wavefunctions 



The quantum field operator hij obeys canonical commutation relations which can be 
derived from the quadratic part of the gravitational action. We have already specified the 



normalization of the creation/annihilation operators (|3.2| ) and of the spatial part of the 



mode functions (|3.3|) . The commutation relations for hij then determines the normalization 
of the time part ipLim{t) of the graviton wavefunctions. The details of this procedure are 
given in and yield a normalization condition 

i^LiUt)jfiimit) - rLimit)j^i^Lim{t) = ?,2n,hG a'^t) . (3.7) 

(Note that the normalization condition given in equation (3.3) of Ford and Parker |^ contains 
a minor typo |T0|.) 



D. Choice Of An Initial (Vacuum) State 

If one defines a Fock vacuum state by the property that it is annihilated by all of the 
operators cum and dumi then the choice of vacuum state is really determined by the choice 
of the mode functions ipLim{t)- For the reasons given in |Q we choose these mode functions 
to be those which correspond to the unique de Sitter invariant vacuum state |0) during the 
initial inflationary stage whose two-point function has Hadamard form. 



E. Wavefunction During The De Sitter Phase 



As shown in equation (2.18) of Ford and Parker 
minimally-coupled massless scalar wave equation 



the mode functions obey the 



'^^ + ^ '^'^ + L{L + 2) 
dt'^ a dt dt 



(3.8) 



There is a slight subtlety: it is impossible to define a de Sitter invariant Fock vacuum state for 
the minimally-coupled massless scalar field |1T[ . However it was shown by Allen and Folacci 
|l^ that the difficulty only arises for the L = mode. In the case of the gravitational field 
operator, the L = and the L = 1 modes are both absent; they correspond to "monopole" 
and "dipole" dynamical degrees of freedom which are not present in the spin-two case. Hence 
in the case of the gravitational field, it is possible to define the desired de Sitter invariant 
vacuum state. The corresponding normalized wavefunction during the de Sitter phase is 



Limit) = ^plit) 



LiL + l){L + 2)\^^ ' adtj 



for t < t, 



end • 



(3.9) 



We note in passing that this time-dependent part of the wavefunction depends only upon 
L and not upon / and m. (This guarantees that the vacuum state will be invariant under 
all rotations of the three-sphere t = constant, which is the subgroup 5*0(4) of the de Sitter 
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group 5*0(1,4)0. However the invariance of the state under the de Sitter group S'0(l,4)o 
is not obvious from inspection. Here the subscript on S'0( 1,4) denotes the part of the 
group connected to the identity.) For this reason, from this point on we drop the indices 
l,m from the time-dependent part of the wavefunction, denoting ipum by ipL- Because 
the wave equation ( p.8| ) is a second-order ODE, the solution ipLit) during the de Sitter 
phase completely determines the solution at all later times. The solution at later times is 
conveniently written in terms of Bogoliubov coefficients. 



F. Wavefunction During The Radiation-Dominated Phase 

The epoch that follows the de Sitter epoch is the radiation-dominated phase. One may 
write the solution to the wave equation during this phase as 

Mt)=aT'rL'{t)+(3l^''rr'{t) for tend<t<teq. (3.10) 

Here, the positive frequency mode during the radiation epoch is defined by 



= -VTT^^"^^"'''^* tend<t<V (3.11) 

a V (L + 1) 

The Bogoliubov coefficients are determined by a condition which follows from the wave 
equation ( p.8|) : both i'lit) and its time derivative must be continuous at all times, and in 
particular at t = tend- One obtains 

ar = {L{L + 2))-'/' (z{L + 1) + ./q^I - ^^^) (3.12) 



prad ^ ^ 1)-\L{L + 2))-l/2gg-2i(L+lK„d_ 

Here Q is the constant defined by 



_ Qoil + Zend)' .3 . 3^ 

^-(fi0-l)(l + Ze,)- ^^-'^^ 

We stress once again that the solution ipL{t) during the de Sitter phase completely determines 
the solution at all later times. In other words the choice of a "positive-frequency" mode 
function during the radiation phase is unimportant. Had we picked a different solution to 
the wave equation ( |3.8| ) to call "positive frequency" then a™'^ and Z?™'^ would have changed 
in such a way as to keep the mode function ipiit) given in ( |3.10D unchanged. In similar 
fashion, the solution of the wave equation during the radiation phase completely determines 
its solution during the matter-dominated phase. 



G. Wavefunction During The Matter-Dominated Phase 

The wavefunction during the matter-dominated (dust) phase may again be expressed as 
a linear combination of the natural positive-frequency solution and its complex conjugate: 
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i.*mat 



for teq<t < 0. 



(3.14) 



The positive frequency mode functions during the matter epoch are 



167rhG 



w/. ^/ . X/ . , , 2i(L + 1) + -^)e"^(^+^)* for tea<t<0. 
aV (^ + 1)(2L + 1)(2L + 3)V ^ ' a dt ) ^ 



The Bogohubov coefficients ai and /5l are determined (as in the spatially fiat case 
combining the Bogoliubov coefficients for the two different phases. 



OIL (3l 
PI «2 



r) \ rad / o \ mat 

OiL Pl\ ( OIL PL 



PI «i 



PI 0^1 



(3.15) 
i)by 

(3.16) 



As previously, the Bogoliubov coefficients and are determined by matching the 
positive frequency radiation mode function ipT'^it) fo the linear combination a^°-^ip^°-^{t) + 
at time teq- One obtains 

((2L + 1)(2L + ?>))-^/^(-2i{L + 1) + ^rw^ + 



„ mat 



iW 



4(L + 1) 



(3.17) 



^mat _ ^ 1)-'((2L + 1)(2L + 3))-i/2p^e-2*(^+^)*"<'. 



Here the constant W is given by 



^^o(l + 



eqj 



no 



(3.18) 



We are now in a position to evaluate the multipole moments {af) of the angular correlation 
function C(7). 



IV. MULTIPOLE MOMENTS OF C{j) 

Combining the results of the previous section, one can easily obtain a formula for the 
multipole moments of the angular correlation function C(7). One replaces the metric per- 
turbation that appears in the Sachs- Wolfe formula ( |2.17| ) with expansion ( [3.1[ ) of the field 
operator. The resulting operator depends upon an angle Q on the celestial sphere. One then 
takes the expectation value of this operator with an identical operator at a different point 
Q' on the celestial sphere. This yields the correlation function 



C(7) 



(4.1) 



1 rO 
4 Jti: 



/° t/t'EL;^ELT„.'^^L(t)^2(tO^^(W)^^'(|i1)>^n^(fi)>^^w(^')(0 



CLlmc\j'l 



. 



Here 7 is the angle between the points Q and Q' on the celestial sphere. Because the Sachs- 
Wolfe formula (p.l7|) involves the time-derivative of the mode function, we have defined 
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■ipL{t) = dipL{t)/dt, where ifji is the mode function during the matter-dominated epoch, 
given in ( |3.14[ ). 



CUmC^L'l' 



To simphfy this expression, first note that the matrix element ^0 

^LL'^iv^mm' ■ This chminatcs the triple sum Y^uvm'- Because the summand is indepen- 
dent of the summation index m, one may then explicitly carry out the sum over m using 
the addition formula for spherical harmonics, equation (3.62) of reference 

' 2/ + 1 

E YUmi'm'^ = ^PKcos7). (4.2) 

m=—l 

Because the argument of the Legendre function Pi{z) is the cosine of the angle 7 between 
the points on the celestial sphere, this shows explicitly that the correlation function depends 
only upon 7. 



1 °° 2/ -I- 1 /"O rO 

= jY.y: ^pkcos7) / dt / dt'MtWLit')R'Lm)R'L{\t'\) m 

4 47r Jh, Ju^ 



Comparing this to the definition of the multipole moments ( |1 . 1| ) , and noting that the sum- 
mation Y1'l=2Ya=2 is equivalent to the summation Y^^2YlT=h o'^^ immediately obtains a 
simple formula for the multipole moment. 



= 7 E / dt dt'^lj,{t)rL{t')R'M)RW\) = 7 E + (4.4) 



1 °° /-O /-O 1 °° 

L=l " *is " ^ L=l 

The complex quantity l\ is what remains of the integral of the mode function along the 
radial null geodesic path. 



\\tR\{\t\)Ll,T\t) (4.5) 

Note that we have assumed (as is implied in Table D that the surface of last scattering lies 
within the matter-dominated epoch; the positive frequency mode function during the matter 
phase is given by ( |3.15|) . The Bolgoliubov coefficients are given by (|3.16|) 

„ „rad„mat , nrad Q*mat „„J o „rad nmat , nrad„*mat /a a\ 

aL = aL + Pl Pl and Pl = aL Pl + Pl (^l > (4-6) 

where the Bolgoliubov coefficients for the matter and radiation transitions are defined by 
( p. 171 ) and ( p.l2| ). In the next section, we discuss how the multipole moments (af) may be 
rapidly evaluated using numerical techniques. 



V. NUMERICAL EVALUATION OF THE MULTIPOLE MOMENTS 

As discussed at the end of Section 0, it is convenient to define dimensionless quantities 
Ml = pjg^""'' '^^g'^'' {af). Using the previous formulae one may write this in the dimensionless 
form 
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where 

J{ = Jl dti?i(|t|) csc2(t/2 + 5)e-^(^+i)* X (5.2) 
-3i{L + 1) cot(t/2 + fi) - I csc2(t/2 + 5) + 2^^ + 41 + 1] 



Taken together with the definitions of and Pl given (|4.(j|), ( ^.IT] ), and ( ^.121 ), the constants 



Q and defined in ( |3.13|) and ( p.l8|) , and the radial function R^ix) defined in ( p.6|) this is 



a self-contained formula for calculating Mi. 

Before discussing the evaluation of Mi in general, it is worth commenting on two limits. 
The first limit is the Zend oo case, where the amount of infiation is large. In this case, 
it is easy to see that Q and hence and Pl diverge oc (1 + .^end)^- Thus in the limit. 
Ml converges. A second interesting limit is the spatially-fiat one, r2o — 1 — >■ 0^, where the 
density parameter approaches unity from above. In this case, it is easy to see that Q and 
hence ai and (3l diverge as {VIq — 1)^^, and the integral diverges as (i^o ~ 1)^^^^. Once 
again, the limit is well-defined. In addition, in this case, the sum over L can be re-written 
as an integral, recovering the k = spatially-fiat formula given in [H]. 



We evaluated Mi using an fourth-order Runge-Kutta adaptive stepsize integrator |T6| to 
obtain the integral which defines Jl- In cases of interest, one frequently needs to include 
many values of L in the summation. In practice we found that summing over the range 
L E 1,1 + I, ■ ■ ■ , /max with /max = 32 + (5/ + 10)/|)f:is| gave results accurate to a few percent for 
reasonable ranges of the free parameters listed in Table |. Rather than compute the J| one 
at a time, it is more practical to compute them "en masse", determining J/, J/^^, ■ ■ ■ , J/^^^ 
simultaneously. This can be done easily because the associated Legendre functions may be 
computed with a stable upwards recursion relation. 



A. Second-Order Recursion Relations 



The upwards recursion relation for the associated Legendre functions is given in equation 
(8.731.2) of reference 



2(1 + j] 



21+ J + i^^(i-i)+m/2l2;j + 2/ + J + l' 0--2)+'+i/2^ 

together with the boundary conditions (or initial values) given in equation (8.755.1) of 

/ sin 



1-J 



-(i+l/2) 



{z) for j = 2, 3, 



(5.3) 



reference ||17|| : 



5-^+1/2), 
1+1/2 ' 



cosx) 



r(/ + 3/2)V 2 



and F;li-i^^\z) 



(5.4) 



These relations may be used to obtain a recursion relation and initial values for the radial 
functions RK ■ The initial values are 
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^'^^) = \l7WTW§^'"^^'''" """^ RUx) = cosx f^f^RKx) (5.5) 
and the recursion relation is obtained from ( |5.3| ) 

for j = 2,3,--- i?l+,(x)= (5.6) 
4a+j)^(/+i-i) „n.,,p< ,'^^1 _ / 0-i)(2;+i)(;+i-i)(;+i-2) / x 

j(i+j+2)(2/+j+l) ^"^A-n-Z+j-lU; Y j(2Z+j+l)(/+j+l)(/+j+2) ^l+j-2\.X) 

Although this recursion relation does not appear to be stable, our experience has been that 
it accurately determines i?^ for I < L < I + 6000. 



VI. NUMERICAL RESULTS AND CONCLUSIONS 

The numerical results are presented as a series of graphs Figs. |I|, 0, and |^ showing 
the values of Mi = ^J'j""'^'^ K^+i) (^af). For all of these graphs, we have taken Z^nd = 10^^, 
Zgq = 10^ and Z\s = 1300, and varied the density parameter ^o- The graphs also show the 
values of Mi for the spatially-fiat k = case, taken from 0. This case corresponds to the 
critically-bound f2o 1 limit. 

It is clear from the figures that this limit is quickly approached; when Qq = 1.1 the Mi 
are almost indistinguishable from the k = spatially-fiat case. It is not hard to see why. 
The effects of the spatial curvature only appear if the past light cone of the observer, taken 
back to the surface of last scattering, actually "probes" a substantial fraction of the spatial 
three-sphere. If the past light cone fails to do this, then within the past light cone the 
universe is indistinguishable (to good approximation) from a spatially-fiat model. 

The fraction of the three-sphere [S^) within this past light cone is easy to determine. 
The three-volume contained within angle Xmax from the point x = of the unit-radius 



may be obtained by integrating ydetPjj where Pij is the three-metric ( |3.4| ). One obtains 

V"(xmax) = 7r(2xmax " sin2xinax)- The total volume of is ^(vr) = 2n'^. If we assume that 
Z\s is much larger than one, then the fraction of the volume of 5''^ contained within the past 
light cone is approximately 

V{\tis\) _ l^(2arcsin 
V{tt) ~ 2^ 
For Qq near 1, this fraction is well-appoximated by 
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nn-i- 



/(fio) = ^ttFt ^ — ^V^^- (6-1) 



f{no)^^iQo-l?. (6.2) 

Thus when Qq = 1.1 the past light cone only explores about 1/1000 of the spatial volume. 
Even if fio = 2 the fraction of the three-sphere that is observed is only f{ilo = 2) = | — ~ 
0.1955 ■ ■ ■. This is why the multipole moments are not very sensitive to provided it is 
close to unity. 

The extension of this calculation to the case of a spatially open FRW universe appears 
straightforward. However it turns out to be much more difficult than expected, primarily 
because the correct choice of initial state is not the obvious one, and because the final result 
for the multipole moments appears to contain logarithmic (infra-red) divergences at zero 
frequency. The spatially open case will be the subject of a forthcoming paper. 
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FIGURES 



FIG. 1. The normalized multipole moments Mi = pj'g."^'' q i'^f) CBR temperature 

fluctuations are shown as a function of the multipole number /, for a spatially-flat (JIq = 1) ^-nd for 
spatially-closed (JIq > 1) cosmological models. — 1 needs to be fairly large for the effects of the 
spatial curvature to be significant. The models being compared all have cosmological parameters 
defined by the redshifts = 1300, Zcq = 10"^ and Zend = 10^^. 

FIG. 2. The normalized multipole moments Mi = ) shown as a function 

of r^o — 1 for I = 2,5, 10, 20, 30, 50. In all cases, the models being compared have the same 
cosmological parameters as in Fig. |l|. Only when Qq becomes significantly larger than one do the 
multipole moments change significantly from the spatially-flat case. 

FIG. 3. The normalized multipole moments Mi = -^|p^^-^-^— ^(of ) are shown as a function of 
r^o ~ 1 for I = 100, 200, 400. In all cases, the models being compared have the same cosmological 
parameters as in Figs. || and||. Only when Qq becomes significantly larger than one do the multipole 
moments change significantly from the spatially-flat case. 
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TABLES 



TABLE I. List of the free parameters that define tlic cosmological model. 



Parameter Units Range Description 

Hq length"^ Ho> Present-day Hubble expansion rate 

dimensionless > 1 Present-day density parameter 

Zig dimensionless Zis > Redshift at last scattering of CBR 

Zcq dimensionless Z^q > Z\^ Redshift at equal matter /radiation energy density 

Zend dimensionless Zend > -^eq Redshift at end of de Sitter inflation 
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